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Abstract 


An improved ability to predict external propulsive performance 
has been incorporated into the three-dimensional Navier-Stokes code 
PA BSD. The improvements are the ability to account for skin fric- 
tion and external pressure forces. Performance parameters for two axi- 
symmetric supersonic cruise nozzle configurations were calculated to 
test the improved methodology. Internal and external flow-field re- 
gions were computed using a two- equation k-e turbulent viscous- stress 
model. The computed nozzle discharge coefficient matched the experi- 
mental data within 0.5 percent in both level and trend. The calculated 
thrust-minus- drag ratios were within 1 percent of the absolute level of 
experimental data , and the trends of data were predicted accurately. The 
predicted trend of integrated nozzle pressure drag matched the trend of 
the integrated experimental pressure drag over a range of nozzle pressure 


ratios, but absolute drag levels were 

Introduction 

Highly maneuverable aircraft operate over a wide 
range of power settings and Mach numbers that re- 
quire a propulsion system with variable geometry for 
obtaining efficient performance at different flight con- 
ditions. Understanding the effects of various nozzle 
geometries on the internal flow region and the sur- 
rounding boattail-nozzlc region is vital for design- 
ing an efficient afterbody for these aircraft. The 
development and utilization of advanced computa- 
tional methods will play a vital role in developing this 
understanding. Several on-going research activities 
currently exist at the Langley Research Center which 
are directed at establishing an experimental database 
for new nozzle concepts. Subsequent improvements 
to the computational methods are guided by these 
data. 

A nozzle internal performance module (ref. 1) has 
been modified to include external aerodynamic ef- 
fects. The module was incorporated into a Navier- 
Stokes solver PAB3D (ref. 2), which provides the 
flow- field solution. The nozzle performance module 
uses the control volume concept to calculate the body 
forces resulting from the fluid flow (ref. 3). Forces 
and moments are calculated from the integration of 
the momentum fluxes through the control volume 
faces. Skin friction on solid walls is calculated by 
using the gradient of the local velocity in the direc- 
tion normal to the wall and the local viscosity. These 
calculations can be performed at intermediate steps 
throughout the solution procedure to provide an inte- 
grated flow quantity convergence in addition to mon- 
itoring the computational residuals in the flow solver 
algorithm. 


ot accurately predicted. 

This paper evaluates the module by using the 
data from the axisymmetric single-engine test body 
with jet- plume simulation reported in reference 4. 
Computed discharge coefficient, thrust- minus-drag 
ratio (measured axial force nondimensionalized by 
ideal isentropic thrust), and internal and external 
surface static pressures were compared with exper- 
imental data for a fixed free-stream Mach number 
of 0.90 at several nozzle pressure ratio settings for 
two of the nozzle concepts in reference 4. 


Symbols 


A e 

nozzle exit area, cm 2 

A-i.j 

local rotation matrix 

Am 

maximum cross-sectional area of 
model, 182.415 cm 2 

At 

nozzle throat area, cm 2 

AA 

incremental cross-sectional 
area, cm 2 

C D,j3 

boattail pressure-drag coefficient 

Q 

discharge coefficient, ^ 

C F ,i 

aerodynamic ideal thrust coefficient, 

Qoc 

Cp.p 

boattail pressure coefficient 

Dri 

nozzle total drag, Pa 

d e 

nozzle exit diameter, cm 

dm 

model maximum diameter, 
15.240 cm 


dt 


F 

F 

Fi 


i,j,k 


l 

l C 


Id 

M 

M 

n 


NPR 

NPR des 

P 

PtJ 

Q 

R 

R 

T 

T tJ 

u 

Ui 

Wi 

w p 

X 


a 

0 


nozzle geometric throat 
diameter, cm 

rate-of-deformation tensor 


7 ratio of specific heats, 1.4 for air 

6 nozzle divergence angle downstream 

of throat, deg 


gross thrust along body axis, N 

total force vector, N 

ideal isentropic gross thrust along 
body axis, N 

indices, when associated with 
computational grid planes, denoting 
orientation in lateral, stream wise, 
and normal directions, respectively 

axial length of boattail, cm 

axial length of nozzle convergent 
section, cm 

axial length of nozzle divergent 
section, cm 

Mach number 

total moment vector, N • m 

local surface normal vector in global 
system 


9 nozzle convergence angle upstream 

of throat, deg 

p viscosity, m 2 /sec 

local Cartesian coordinate system 
on surface 

p density, kg/m 3 

T{j viscous-stress tensor, Pa 

<t>i local Cartesian skin friction force 

vector, N 

Subscripts: 

oo free-stream conditions 

flux flux and pressure forces (flow- 

through face) 

pres pressure forces (solid face) 

fric viscous forces (solid face) 

ref reference conditions 

local local conditions 


nozzle pressure ratio, ^ 

design nozzle pressure ratio 

static pressure, Pa 

jet total pressure, Pa 

dynamic pressure, Pa 

gas constant (7 = 1.4), 

287.3 J /kg * K 

moment arm vector, cm 

static temperature, K 

jet total temperature, K 

total velocity vector, m/sec 

local Cartesian velocity, m/sec 

ideal mass-flow rate, kg/sec 

actual mass- flow rate, kg/sec 

axial distance, positive down- 
stream, cm 

nondimensional distance from wall 
in turbulent boundary layer 

angle of attack, deg 

nozzle boattail angle, deg 


Abbreviations: 

A/B afterburning 

C-D convergent-divergent 

CFL Courant-Friedrichs-Lewy 

C2 configuration 2 of reference 4 

C4 configuration 4 of reference 4 

2- D two-dimensional 

3- D three-dimensional 

Experimental Model Configuration 

General Model Description 

The experimental results selected for this com- 
putational study are reported by Carson and Lee 
in reference 4. A photograph of the model in the 
Langley 16-Foot Transonic Tunnel and a sketch of 
a representative test model and support system are 
presented in figures 1 and 2. The facility, which is a 
single-return, continuous-flow atmospheric wind tun- 
nel with an octagonal slotted-throat test section, has 
a continuously variable Mach number range from 0.20 
to 1.30. A detailed description of this wind tunnel is 
given in reference 5. 

The models of the axisymmetric, convergent- 
divergent (C-D) nozzles (figs. 3 and 4) had a circular- 
arc throat contour and conical divergent sections. 
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The nozzle geometries simulated a variable-geometry 
axisymmetric nozzle designed for a variable- cycle en- 
gine that was proposed for a supersonic cruise air- 
craft. Experimental data were obtained for five noz- 
zle configurations that represented the geometries for 
different flight conditions and power settings. These 
geometries ranged from a subsonic cruise, dry power 
configuration, which would have a low expansion ra- 
tio and a high nozzle boattail angle, to a super- 
sonic acceleration, maximum afterburning configu- 
ration, which would have a high expansion ratio 
and a low nozzle boattail angle. Tests were con- 
ducted at nozzle pressure ratios (NPR’s) from jet 
off to approximately 10 at free-stream Mach num- 
bers from 0.60 to 1.30. The nozzles were attached 
to an axisymmetric single-engine host body mounted 
on a sting-strut system and had a maximum cross- 
sectional area of 182.415 cm 2 (28.274 in 2 ). The ogive- 
shaped nose had an apex angle of 14° and a radius 
of curvature of 128.37 cm (50.54 in.). The nonmetric 
portion of the host body was 67.31 cm (56.50 in.) 
long. The model was cylindrical from the metric 
break to the nozzle connect station and was 69.85 cm 
(27.50 in.) long. 

The jet was simulated by high-pressure air ex- 
iting the nozzle into the external free stream. An 
air system at the wind tunnel provides a continuous 
flow of clean, dry air at a stagnation temperature 
of nominally 300 K (540°R). The nozzle discharge 
coefficient was determined from experimentally mea- 
sured jet total temperature, jet total pressure, and 
measured mass-flow rate. The thrust ratio was deter- 
mined from the measured balance axial force that was 
nondimcnsionalized by the ideal thrust (determined 
from the measured nozzle mass flow) and corrected 
for internal pressure tare forces. Static pressure ori- 
fices were located on the external boattail and the 
internal nozzle surfaces. 

Test Case Details 

Two nozzle configurations from reference 4 were 
selected as test cases for the computational method, 
and they are described below. 

Configuration 2 . Configuration 2 of reference 4 
is a high expansion ratio, axisymmetric, C-D nozzle. 
The external geometry is a circular- arc and straight- 
line type with a terminal boattail angle of 3.82°. The 
internal geometry has a circular-arc throat and a di- 
vergent section angle of 13.18° with an expansion ra- 
tio of 3.000 that results in a design NPR of 21.23. 
This geometry was chosen because the experimental 
data generally indicate no external flow separation 
at Mso = 0.90 for the NPR settings tested. In spite 


of the absence of data around the design NPR at 
Mryj = 0.90 because of limitations in the attainable 
mass- flow rate of the experimental apparatus, it re- 
mains an interesting case for evaluation of a compu- 
tational method. The internal flow does separate at 
lower NPR settings because of the high expansion 
ratio. 

Configuration 4* Configuration 4 of reference 4 
is a low expansion ratio, axisymmetric, C-D nozzle. 
The external geometry is a circular-arc and straight- 
line type with a terminal boattail angle of 8.28°. The 
internal geometry has a circular- arc throat and a 
straight-line divergent section angle of 4.78° with an 
expansion ratio of 1.500 and a design NPR of 6.23. 
This geometry was chosen because the experimental 
data were obtained around the design NPR and 
nozzle internal flow appears to remain attached at the 
lower NPR conditions because of the low expansion 
ratio. 

Computational Procedure 

Flow-Field Calculation 

A three-dimensional Navier-Stokes code PAB3D 
was developed to predict the effects of the jet exhaust 
plume on nozzle- afterbody configurations. The thin- 
layer Navier-Stokes formulation (ref. 6) was modified 
to simulate jet mixing problems (ref. 2). The code 
allows for the discretization of the flow-field domain 
into multiblock grids and can utilize several numeri- 
cal schemes to solve the governing equations and tur- 
bulence models as discussed in reference 7. The Roe 
upwind scheme with third-order accuracy is used to 
evaluate the explicit part of the governing numerical 
equations for relaxation in the streamwise direction, 
typically along the i index. The van-Leer scheme, 
with its faster convergence rate, is used to construct 
the implicit operator in the cross plane; typically this 
plane is a grid along the j and k indices. (Normally 
this code utilizes the Roe scheme to sweep streamwise 
through the computational domain and the van-Leer 
scheme for the solution of the cross plane.) For the 
present study, the streamwise grid plane is oriented 
along the j and k indices to take advantage of the im- 
plicit scheme for efficiency and accuracy. Solutions 
were developed using two versions of a two-equation 
k-e turbulent viscous-stress model. A low Reynolds 
number k-c model was used for the flow near solid 
surfaces while a high Reynolds number k-e model was 
used in the regions of free-shear flow. 

A user-written control file determines the com- 
munication between blocks and the type of bound- 
ary condition to be used at each face. Different grid 
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topologies for neighboring blocks and mixed bound- 
ary conditions on a block face arc permitted, with 
some restrictions on grid matching at block bound- 
aries. The code permits different numerical schemes 
to be selectively applied to each block. 

Grid Definition 

The block arrangement and dimensions of the 
grid are shown in figure 5(a). A sectional view of 
the overall grid in the j and k planes is shown in 
figure 5(b), and a sectional view of the two blocks 
defining the internal flow path is shown in figure 5(c). 
A 2-D wedge grid was used, and the flow was assumed 
to be axisymmetric. A single-cell- width cylindrical 
wedge grid was generated with a wedge angle of 5.00°, 
which represents V 72 of an axisymmetric geometry. 

The nozzle flow field was modeled by using five 
computational blocks. Two blocks represented the 
nozzle internal grid, and an additional two blocks 
represented the external grid region surrounding the 
nozzle. The fifth block modeled the external re- 
gion downstream of the nozzle exit. The dimen- 
sions of the nozzle internal blocks were 2 x 6 x 49 
and 2 x 78 x 49. The nozzle internal boundary-layer 
grid region contained approximately 24 points with 
the first grid point specified by — 2.5. The ex- 
ternal blocks surrounding the nozzle grid were dimen- 
sioned 2 x 6 x 87 and 2 x 64 x 87. Approximately 
40 streamwise grid points were used in the region of 
the nozzle boattail. The external downstream region, 
modeled by block 5, was dimensioned 2 x 51 x 135. 
The interfaces matched point for point between the 
blocks upstream and downstream of the nozzle exit 
plane. The grid was extended to a distance of 15 noz- 
zle exit radii away from the outer boundary. The 
outflow boundary was 5 exit radii downstream of the 
nozzle exit plane. 

Boundary Conditions 

The PAB3D code allows different boundary con- 
ditions to be applied to a given block face. Solid walls 
are treated as no-slip adiabatic surfaces. The inflow 
boundary conditions used for the internal nozzle flow 
path are the total pressure pf j and the total temper- 
ature T t j . This particular inflow boundary condition 
assumes a fluid flow angle normal to the inflow face. 
The operating NPR of the nozzle and the free-stream 
static pressure c determines the jet total pressure 
Ptj = (NPR)(Pbo)- 

A jet total temperature of 300 K was used for all 
the calculations, and it was the nominal temperature 
set during the experimental investigation. The in- 
flow Mach number, flow angle, total pressure, and to- 
tal temperature were fixed along the external inflow 


boundary. A boundary condition for Riemann invari- 
ants (ref. 2) along the characteristics was specified for 
the lateral free-stream boundary of the flow domain. 
An extrapolation boundary condition was applied on 
the downstream outflow face where both the exter- 
nal free stream and the nozzle plume stream exit the 
computational domain. The wedge- angle boundary 
condition determines the angle of the wedge lateral 
faces and the proper boundary values to apply for 
satisfying the symmetry assumptions. Flow proper- 
ties were transferred between blocks through cell-by- 
cell conservation of the mass and momentum fluxes 
across the block interface with third-order continu- 
ity. However, some restrictions on the multiblock 
grid were expected by the code, such as an integer 
multiple to one correspondence between cells at the 
block interface and matched cell sizes in directions 
normal to the block interface. 

Performance Calculation 

Nozzle performance is obtained through the ap- 
plication of the momentum theorem to a control 
volume surrounding the nozzle (ref. 3). Cheatham, 
Walker, and Gridley calculated both 2-D and 3-D in- 
viscid nozzle performance in reference 8, and Carlson 
calculated 3-D viscous nozzle performance in refer- 
ence 1 by using the control volume method. Inte- 
gration of the flow quantities is typically performed 
across the nozzle exit. The method integrates the 
mass and momentum fluxes and the pressure forces 
over each cell by using equations (1) and (2) for flow- 
through sections of the control volume: 

W P = J2 AA C 1 ) 

Fflux = ^2 [pu(u - n) + (p - Poo)n] A.4 (2) 

where A A is the incremental cross-sectional area of 
the cell face and R is the moment arm vector from 
the reference center to the cell face. 

Ideal mass-flow rate and thrust are determined 
from the isentropic flow equations (3) and (4), re- 
spectively (ref. 9), and they are used to normalize 
the calculated mass-flow rate and thrust for compar- 
isons with the experimental data: 
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( 11 ) 


Skin friction and pressure forces are calculated for 
solid wall sections of the control volume. The solid 
surface pressure force is calculated by extrapolating 
the cell centered static pressure to the surface wall 
and multiplying this pressure by the cell face area: 


F — Fflux T F pres + Ff r j c 

The performance package is incorporated into 
PAB3D to permit monitoring of various performance 
parameters as the solution convergences. 


Fpres — ^ ^ [( P Poo)fi] AT 


Discussion of Results 


The viscous-stress tensor used for determining the 
skin friction force is calculated by using only the ve- 
locity derivatives normal to the surface. The velocity 
gradients are determined by a two-point difference. 
The first velocity is a zero magnitude vector posi- 
tioned on the surface. The second velocity is the 
velocity at the cell center as sketched in figure 6. 
Equations (6a) and (6b) are two of the nonzero com- 
ponents of the rate-of-deformation tensor calculated 
(ref. 10): 


e 13 


623 


dU\ 

dU2 

9^3 


(6a) 

(6b) 


where U{ — AijUj. The remaining velocity deriva- 
tives are assumed to be zero. The set of local velocity 
components Ui are in a local surface Cartesian coor- 
dinate system & (where £3 is the direction normal to 
the surface) , which is determined from an orthogonal 
transformation of the global velocity components by 
the rotation matrix A r j . In general e^j = where 
i ± 3- 


The local shear-stress tensor where i / j, 
was the local viscosity multiplied by the velocity gra- 
dients eij , where i 7^ j. The viscosity was deter- 
mined from Sutherland’s formula (ref. 11) by using 
the static temperature at the local cell center: 


AMocal — /Tef 


(T ref + 119.33) 


1.50 


(Tlocal + 1 10.33) 


( 1 ?) (7) 


Experimentally measured external-surface static- 
pressure coefficient distributions and the internal- 
surface pressure ratio distributions are compared 
with calculations in figures 7 to 12. Various measured 
force parameters are compared with calculations in 
figures 13 to 17. Specific comparisons are presented 
in the following figures: 


Figure 


External and internal surface distributions 

of configuration 2, NPR = 4, Mqq = 0.90 . . . 7 
External and internal pressure distributions 

of configuration 2, NPR = 5, Moo = 0.90 . . . 8 
External and internal pressure distributions 

of configuration 2, NPR = 6, Moo — 0.90 ... 9 
External and internal pressure distributions 

of configuration 4, NPR = 5 ; Moo = 0.90 . . 10 

External and internal pressure distributions 

of configuration 4, NPR = 6, A/00 = 0.90 . . 11 

External and internal pressure distributions 

of configuration 4, NPR = 7, M 00 = 0.90 . . 12 

Nozzle discharge coefficients 13 

Variation in nozzle ideal thrust coefficient 

with NPR 14 

Variation in aeropropulsive performance with 

NPR 15 

Variation of increment in aeropropulsive 

performance with nozzle configuration ... 16 

Variation in integrated nozzle pressure drag 

with NPR 17 


Discussion 


Tij = Mocal e ij (* # j) (8) 

The local skin friction force components 0, were 

4>i = ( r ij n i) AA (9) 

where n* is the ith component of the unit normal 
vector n and 4> t were transformed back to a force 
vector expressed in global coordinates F f r j c . The 
solid surface forces are then added to Ffl ux for the 
total volume forces: 

(Ffri c )j - ^Aj. (10) 


Solutions were obtained for two nozzle configu- 
rations by using the PAB3D flow solver. The solu- 
tions were computed by using a single-cell-width 2-D 
mesh oriented so that the entire flow field was solved 
implicitly with each iteration. All solutions were 
developed by using a two-equation, k-e turbulent 
viscous-stress model. Discharge coefficient, ideal 
thrust coefficient, thrust- minus-drag ratio, and noz- 
zle drag coefficient were calculated for different noz- 
zle operating conditions, which exhausted into a 
Moo — 9*99 free-stream airflow at an angle of attack 
of 0°. The calculations are summarized as follows: 
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NPR 

Parameter 

4 

5 

6 

7 

Configuration 2 

Iteration ...... 

8000 

6800 

2000 

2500 

Computer time, hr 

3.3 

2.8 

0.8 

1.0 

C d 

0.963 

0.963 

0.963 

0.963 

(■ F-Dn)/Fi .... 

0.789 

0.787 

0.849 

0.887 

Configuration 4 

Iteration 


2500 

2500 

1500 

Computer time, hr 


1.0 

1.0 

0.6 

c d 


0.966 

0.966 

0.966 

(F-D n )/Fi .... 


0.969 

0.979 

0.984 


The solutions started with a CFL number of 1. 
The CFL number for subsequent iterations during 
a run were adjusted automatically within the code 
for convergence acceleration. Flow solution residuals, 
discharge coefficient, thrust ratio, and nozzle drag 
were obtained at regular intervals during the solution 
development. Typical solution convergence criteria 
are less than 0.05 percent change in discharge coef- 
ficient, less than 0.1 percent change in thrust ratio, 
and less than 0.001 change in nozzle drag coefficient 
over the previous 1000 iterations. Solutions with flow 
separation typically required an additional 3000 to 
4000 iterations beyond the average 1500 to 3500 iter- 
ations to establish an acceptable level of performance 
convergence. Situations of flow instability because of 
either real flow physics or numerical modeling diffi- 
culties occasionally caused some performance param- 
eters to oscillate about a median value. The solutions 
for this study met the general convergence criteria 
with little or no occurrence of oscillations in perfor- 
mance numbers. 

Comparisons of predicted internal and exter- 
nal static pressure distributions with experimental 
data are made for configuration 2 at NPR = 4, 5, 
and 6. Predicted static pressure distributions arc 
compared with experimental data for configuration 4 
at NPR = 5, 6, and 7. Predicted performance pa- 
rameters are compared with experimental data. 

Static-Pressure Comparisons 

Computed external and internal static-pressure 
distributions are compared with experimental data 
for configuration 2 in figures 7 to 9. In general, the 
external static-pressure coefficients are matched well 
for all three NPR conditions tested (figs. 7(a), 8(a), 
and 9(a)). The corner flow expansion near = 0 
is slightly overpredicted and the trailing-edge flow 


expansion is slightly low, which is possibly a result of 
grid density in those locations. The flow separation 
over the internal divergent section near = 0.7 is 

®77l 

fairly closely matched (fig. 7(b)) at NPR = 4. Neither 
the apparent internal flow separation at NPR = 5 
(fig. 8(b)) at -r- = 0.92 nor the trailing-edge flow 
compression at NPR = 6 (fig. 9(b)) at = 0.995 
was predicted. Extremely fine grid in the vicinity of 
the nozzle trailing edge in the external and internal 
regions would probably be required to resolve such 
detailed local flow structure. The slight mismatch of 
the computed pressures and the experimental data 
upstream of — 0.20 is probably because of a 
geometry documentation difficulty in reference 4. 
Analytic results performed in reference 4 and the 
present method showed very similar trends in the 
internal static- pressure distributions in the region 
upstream of the throat. 

Computed external and internal static-pressure 
distributions for configuration 4 are compared with 
experimental data in figures 10 to 12. The predicted 
external pressure-coefficient distributions match well 
with experimental data except for the magnitude of 
the flow expansion around the circular-arc clqsure of 
the boattail ^0 < ^ < 0.2^) in figures 10(a), 11(a), 

and 12(a). The overprediction of this flow accelera- 
tion is, again, possibly because of grid density in that 
region of the model. The internal flow for this nozzle 
remained fully attached for the range of NPR ex- 
amined. The internal static-pressure ratio distribu- 
tion downstream of the throat typically changed little 
for choked nozzle flow that filled the internal nozzle 
volume completely. The solutions matched the ex- 
perimental internal pressure distributions closely for 
the three NPR computed and shown in figures 10(b), 
11(b), and 12(b) (NPR = 5,6, and 7, respectively). 

Performance Comparisons 

Figures 13 to 16 show comparisons of com- 
puted performance parameters with experimental 
data from reference 4. All data were at a free-stream 
Mach number of 0.90 and a free-stream angle of at- 
tack of 0°. Additional experimental data were plot- 
ted to show general trends that surrounded the con- 
ditions chosen for the analytic study. 

Discharge coefficient A comparison of pre- 
dicted discharge coefficients with experimentally 
determined discharge coefficients is shown in fig- 
ure 13. The experimental data, shown by the circle 
and square symbols, were obtained at static condi- 
tions. Theoretically, external flow has no effect on 
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discharge coefficients once the internal flow is choked 
(NPR > 1.89). Thus, a comparison of computed dis- 
charge coefficients at Moo = 0.90 with experimental 
discharge coefficients at Moo — 0.0 is valid. The pre- 
dicted discharge coefficients for configuration 2 are 
within 0.5 percent of the experimental data. The pre- 
dicted discharge coefficients for configuration 4 are 
within 0.1 percent. 

Ideal thrust coefficient . Predicted ideal thrust 
coefficients for both configurations are compared 
with experimental data in figure 14. The ideal 
thrust coefficient is the ideal nozzle thrust non- 
dimensionalized by using the free-stream dynamic 
pressure goo, a reference area A m , and the physical 
mass flow w p . This coefficient is used as a means 
of converting aerodynamic coefficients normalized by 
Qoc A m to propulsive- force ratios based on Fj. For 
these configurations, the reference area is the maxi- 
mum cross-sectional area of the model. The ability to 
predict this quantity should be similar to the ability 
to predict the discharge coefficient. In general, the 
level and slope of the computed ideal thrust coeffi- 
cient and the experimental data are closely matched, 
which indicates that the trends in the mass flow with 
NPR were predicted properly. 

Aeropropulsive performance . Figure 15 

shows a comparison of the predicted aeropropulsive 
performance (thrust-minus-drag ratio) with the ex- 
perimental data for both configurations for 
Moo = 0.90. The flow that exhausted from the 
high expansion ratio nozzle (configuration 2) is highly 
overexpanded and exhibits some internal flow sepa- 
ration at NPR = 4, which causes the discontinuity in 
the trend between NPR = 4 and 5 (square symbols, 
fig. 15). The method predicted the change in perfor- 
mance as the region of nozzle internal flow separation 
diminished. The predicted thrust-minus-drag ratio 
was within 0.010 of the experimental data. Consid- 
erably higher performance was attained by configu- 
ration 4 as the design NPR was bracketed by the 
NPR range from 5 to 7. The level and trend of the 
predicted thrust-minus-drag ratio for configuration 4 
was within 0.004 of the experimental data in thrust 
ratio, when the circle symbols are compared with the 
solid line. 

The difference in aeropropulsive performance be- 
tween configurations 2 and 4 was predicted within 
0.005 in thrust ratio for the three NPR settings 
shown in figure 16. Limitations in calculating the 
internal flow separation of the turbulence model for 
configuration 2 likely contributed to the disagree- 
ment with the experimental data at NPR — 5. The 


effect on performance of any flow separation, which 
may be present at NPR = 6 and 7, appeared to be 
less significant than at NPR = 5. The comparison 
of calculated with experimentally determined incre- 
ments in performance was within 0.002 in thrust ratio 
for NPR = 6 and 7. 

Nozzle pressure- drag coefficient . The inte- 
grated nozzle pressure-drag coefficient is shown in 
figure 17. The trend in pressure-drag coefficients 
with NPR is closely predicted for both configura- 
tions, although the absolute level of drag is not 
well matched. Predicted pressure-drag coefficients 
for configuration 2 are on the average 0.006 below 
the experimentally determined pressure-drag coeffi- 
cients. The opposite condition occurred for config- 
uration 4 where the predicted pressure drag is typi- 
cally 0.004 above the experimentally determined drag 
coefficients for the three NPR conditions examined. 
The lower than predicted drag of configuration 2 was 
possibly caused by the underprediction of the flow 
expansion near the nozzle trailing edge. Similarly, 
the higher than predicted drag of configuration 4 was 
possibly caused by the overprediction of the expan- 
sion on the nozzle boat tail shoulder. Accurate pres- 
sure drag numbers are sometimes difficult to obtain 
from wind-tunnel models through the pressure-area 
integration method. The ability to predict trends, 
however, is very useful. The prediction of accurate 
levels of performance may occur with improvements 
in both testing and computational techniques. 

Conclusions 

A nozzle internal performance module, which is 
part of the three-dimensional Navier-Stokes PAB3D 
method, has been modified to include external aero- 
dynamic flow effects. The flow quantities calcu- 
lated were external and internal surface static pres- 
sure, discharge coefficient, ideal thrust coefficient, 
thrust-minus-drag ratio, and pressure-drag coeffi- 
cient. Comparisons of the results of the Navier- 
Stokes PAB3D method were made with a selection of 
experimental data from a model of an axisymmetric 
supersonic cruise nozzle investigation performed in 
the Langley 16-Foot Transonic Tunnel. Calculations 
were performed for a single transonic Mach number 
at several nozzle pressure ratio settings. Solutions for 
operating conditions close to design and below design 
nozzle pressure ratios were obtained. The results are 
summarized as follows: 

1. Internal and external surface static-pressure dis- 
tributions were closely predicted for the selected 

configurations. In particular, the location of the 
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internal separation of flow from the divergent sec- 
tion of the high expansion ratio configuration at 
a nozzle pressure ratio of 4 was predicted closely. 

2. Nozzle discharge coefficient was predicted to 
within 0.5 percent of experimental data for nozzle 
pressure ratios at design and overexpanded nozzle 
flow conditions. The trends and the levels of ideal 
thrust coefficients matched the experimental data 
well. 

3. Thrust-minus-drag ratios were predicted to with- 
in 0.010 in thrust ratio. The predicted increment 
in thrust-minus-drag ratio between the two noz- 
zles for nozzle pressure ratios of 5, 6, and 7 was 
within 0.005 of the experimental data. 

4. Integrated nozzle pressure-drag coefficient trends 
with nozzle pressure ratio matched fairly closely; 
however, absolute levels of drag coefficient were 
not well predicted. 

NASA Langley Research Center 

Hampton, VA 23681-0001 

November 12, 1993 
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Figure 1. Photograph showing installation of a typical nozzle in test section of Langley 16- Foot Transonic Tunnel. 
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Figure 2. General arrangement of model and support system. AH linear dimensions are in centimeters. (From ref. 4.) 
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Figure 3. Photograph showing five nozzle configurations tested (ref. 4). 
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Figure 4. Geometric details of test nozzle configurations. Absolute linear dimensions are in centimeters. 
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(a) Blocking 
Figure 5. 
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(b) Internal surface. 


Figure 7. Comparison of calculation with experimental data for surface-pressure distributions of nozzle 
configuration 2 with NPR = 4 and Moo = 0.90. 






P/Ptj 



Figure 8. Comparison of calculation with experimental data for surface-static pressure distributions of nozzle 
configuration 2 with NPR = 5 and Mqc = 0.90. 
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(b) Internal surface. 


Figure 9. Comparison of calculation with experimental data for surface-pressure distributions of nozzle 
configuration 2 with NPR = 6 and = 0.90. 
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(b) Internal surface. 


Figure 11. Comparison of calculation with experimental data for surface-pressure distributions of nozzle 
configuration 4 with NPR = 6 and M x — 0.90. 
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Figure 12. Comparison of calculation with experimental data for surface- pressure distributions of nozzle 
configuration 4 with NPR — 7 and M 0 0 — 0.90. 
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Experiment PAB3D Configuration 
M = 0.90 (ref. 4) M = 0 (ref. 2) 



NPR 

Figure 13. Comparison of calculation with experimental data for nozzle discharge coefficients. 


Experiment PAB3D Configuration 

(ref. 4) (ref. 2) 
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Figure 14. Comparison of calculation with experimental data for variation in nozzle ideal thrust coefficient 
with NPR with M 0 0 = 0.90 and a = 0°. 
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Experiment PAB3D Configuration 

(ref. 4) (ref. 2) 


□ 2 



Figure 15. Comparison of calculation with experimental data for variation in aeropropulsive performance with 
NPR with Moo = 0.90 and a = 0°. 



NPR 

Figure 16. Comparison of calculation with experimental data for variation of increment in aeropropulsive 
performance with nozzle configuration. 
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Experiment PAB3D Configuration 

(ref. 4) (ref. 2) 

O 4 

□ 2 
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Figure 17. Comparison of calculation with experimental data for variation of integrated nozzle pressure drag 
with NPR with M, ^ = 0.90 and a = 0°. 
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